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ABSTRACT 



We consider the precession of isolated neutron stars in which supcrfiuid is not pinned to 
the stellar crust perfectly. In the case of perfect pinning, Shaham (1977) showed that there are 
no slowly oscillatory long-lived modes. When the assumption of perfect pinning is relaxed, 
new modes are found that can be long-lived, but are expected to be damped rather than 
oscillatory, unless the drag force on moving supcrfiuid vortex lines has a substantial component 
perpendicular to the direction of relative motion. The response of a neutron star to external 
torques, such as the spindown torque, is also treated. We find that when computing the response 
of a star to perturbations, assuming perfect coupling of supcrfiuid to normal matter from the 
start can miss some effects. 
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1. INTRODUCTION 

Radiopulsars can be exceptionally stable clocks. The predictability of pulse arrival times has made 
possible precision tests of General Relativity (e.g. Taylor et al. 1992) and the discovery of the first 
extrasolar planetary systems (Wolszczan and Frail 1992). However, imperfections in the rotation of most 
pulsars have been monitored for some time, most notably glitches (e.g. Boynton et al. 1969, Radhakrishnan 
and Manchester 1969, Backus et al. 1982, Downs 1982, Demianski and Proszynski 1983, Manchester et al. 
1983, Lyne 1987, Lyne and Pritchard 1987, Cordes, Downs and Krause-Polstorff 1988, McKenna and Lyne 
1990, Hamilton et al. 1989, McCulloch et al. 1990, Flanagan 1990, 1993, Lyne, Smith and Pritchard 1992, 
Shcmar and Lyne 1996) and timing noise (e.g. Boynton et al. 1972, Groth 1975, Cordes and Helfand 1980, 
Cordes and Downs 1985, D'Alessandro et al. 1995). These deviations from stable spin convey information 
about the internal structure and dynamics of neutron stars. For example, the long-time healing of the spin 
frequencies and spindown rates of glitching pulsars may be explained if the convulsions are due to sudden 
unpinning of superfluid vortex lines that were glued to nuclei in the pulsar's crust before the glitch, migrate 
as a consequence of the glitch, and repin to other nuclei in its aftermath (e.g. Anderson and Itoh 1975, 
Alpar et al. 1981, 1984a,b, 1993, Link, Epstein and Baym 1993). 

A small number of neutron stars also exhibit long term cyclical but not precisely oscillatory variations 
in their spin. A particularly well-known case is the Crab Pulsar, whose phase residuals (after careful fitting 
that accounts for spindown and glitches) vary systematically with a peak-to-peak range of order ±10 ms, 
and a characteristic cycle duration of about 20 months. (Lyne, Pritchard and Smith 1988). After its 
Christmas 1988 glitch, the Vela Pulsar showed - after accounting for exponential recovery from the glitch 
damped oscillatory phase residuals with a period of order 25 days; evidence for oscillations in the frequency 
derivative of the pulsar both before and after the glitch with a period of "a few tens of days" was also 
reported (McCulloch et al. 1990). Evidence for long term variations (correlation times ~ 100 days) in the 
pulse shape of the Vela Pulsar has also been found in data spanning approximately four years (Blaskiewicz 
1992, Cordes 1993). A principal component analysis of the pulse shape of PSR 1642-03 (Blaskiewicz 1992; 
Blaskiewicz and Cordes 1997, in preparation) yields evidence for cyclical pulse shape variations with a 
period of about 1000 days; long-term variations on a similar characteristic timescale are also seen in the 
timing residuals for this pulsar (Cordes 1993). Finally, although Her X-l is an accreting X-ray pulsar, not 
an isolated radiopulsar, it has a well-known 35 day cycle on which it appears and disappears; observed 
variations in pulse shape over the cycle suggest that it is related to periodic variations in the rotation of the 
neutron star (e.g. Triimper et al. 1986, Alpar and Ogclman 1987). 
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Soon after the discovery of radiopulsars, it was suggested that long term variations in their spin could 
result from free precession (e.g. Davis and Goldstein 1970, Goldreich 1970, Ruderman 1970, Brecher 1972, 
Pines et al. 1973, Pines and Shaham 1972, 1974, Lamb et al. 1975, Jones 1976), and that pulsar arcanac 
such as drifting subpulses might be related to preccssional effects. If a neutron star were a rigid or semirigid 
solid (see Pines and Shaham 1972, 1974), it would precess with a period of order P/e, where P is the spin 
period of the star and e is the fractional difference between its principal moments of inertia; e < 10~ 7 would 
imply precession periods >P(s) years, where P(s) is the spin period in seconds. However, once it became 
apparent that superfluid vortex lines pin to the crust of a neutron star, Shaham (1977) demonstrated that 
slow, persistent precession is impossible. Pinned superfluid alters the effective asymmetry of the star to 
Ip/I, where I p is the moment of inertia of the pinned superfluid and I the moment of inertia of the star (or 
stellar crust, depending on various coupling parameters); since Ip/I ~ 10~ 2 or even larger, precession is very 
fast. Moreover, Shaham (1977) demonstrated that precession would decay rapidly for some estimates of the 
coupling timescale between crust and core (e.g. Alpar and Sauls 1988; Sedrakian and Sedrakian 1995). 

Perhaps because of Shaham's (1977) pessimistic conclusions, there has been relatively little theoretical 
work on the long term variability of pulsar spins. Some argue that vortex line pinning does not occur, 
making slow precession possible (e.g. Jones 1988 for the Crab Pulsar), but it may be hard to support the 
viewpoint that pinning is completely absent in face of physical arguments to the contrary (e.g. Alpar et al. 
1984a, Epstein and Baym 1988, Link and Epstein 1991). Notwithstanding the pessimism of theorists, the 
data demand an explanation. Ruderman (1997) has mentioned long term variations in pulsar spin rates as 
one of the outstanding unresolved problems of neutron star physics. 

This paper is the first of two that attack this problem. Our purpose in this paper is three fold: First, 
we revisit the arguments put forward by Shaham (1977) with an eye toward identifying possible loopholes. 
Although it will become apparent that there are several different possibilities, we concentrate here on 
Shaham's assumption that superfluid pins perfectly to crustal nuclei. (Some of the other loopholes will 



be considered in a subsequent paper.) We develop the formalism for doing so in Section 2.2, and solve 
the equations governing the spin dynamics of a neutron star in succeeding sections to varying degrees of 
complexity and realism. As expected, we do find modes in addition to those found by Shaham (1977), but 
we also argue that none of these modes is likely to be a long period, slowly damped oscillation. Second, we 
examine what can happen in a multicomponent star, in which some regions contain pinned superfluid, and 
others unpinned superfluid. One might think that if some parts of a star are capable of slow oscillations 
of the spin - either precession or long period fluid modes, such as Tkachenko modes - then there could 
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be an observable signature of these modes in the detectable pulsar spin rate. However, we demonstrate 
that this situation is highly unlikely, for even if such regions do exist inside actual pulsars, persistent long 
period oscillations in those domains are only possible if the coupling to the crust, where superfluid is pinned 
effectively, is very weak; but under such conditions, the crust is almost unaffected by the slow oscillations, 
which hardly manifest themselves in the crustal spin rate. Third, we begin a general examination of the 
effects of external torques - such as the spindown torque - on the spin dynamics. For this purpose, we 
derive explicit expressions for the response of the various components of a neutron star to rather general 
time dependent torques. Our treatment of this problem shows that the limit of perfect coupling must be 
taken carefully when the response to external torques is needed, because the additional modes that appear 
when pinning is imperfect contribute to the response, and cannot be ignored. 

As will become apparent in the succeeding sections, we do not believe that the long term cyclic 
variability detected in the spins of some pulsars can be accounted for by free precession, and that it is not 
likely to be due to forced precession either. However, we do believe that this paper begins to elucidate the 
complexity of the behavior of neutron star spin, and clarifies the conditions that must be met for precession 
to occur, even if those conditions are not likely to be realized. 



2. OVERVIEW 

2.1. Pinned Superfluid Suppresses Precession (Shaham 1977) 

Shaham (1977) showed that pinned crustal superfluid dramatically alters the physics of precession. Let 
us review his argument briefly. Consider a three component neutron star that consists of: (i) a rigid crust 
rotating at angular velocity i~2 cr ; (ii) pinned crustal superfluid, whose angular momentum L p is independent 
of time in the frame rotating with the crust; and (iii) a core (super)fluid rotating at angular velocity fi c . 
As seen in the inertial frame, 



dfl c d(I cr ■ n cr ) 

i c — + — - — + n cr xL p = o, (i) 



if there are no external torques, where I c is the moment of inertia of the core fluid and I cr is the moment of 
inertia tensor of the crust. We assume that the moment of inertia tensor of the core fluid is always of the 
form (S is the unit tensor) 

I c = (I c AI C )S + A/ c ( ™ c nc-5 \ ^ (2) 
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so that its angular momentum L c = I c ft c = I c ft c . This amounts to assuming that the core fluid adjusts 
its shape instantaneously to an oblate spheroid flattened along its direction of rotation. We shall discuss 
this assumption more fully in a subsequent publication. 

Introducing a dissipative torque that seeks to enforce corotation between the crust and the core, we get 
the coupled equations 



^(Icr ' ^cr) 



dt 



+ fl cr X L r 



-K(n cr - n c ) 

«c), 



(3) 
(4) 



where K is a constant. These equations have a rich set of fixed points (where time derivatives vanish) 
depending on the orientation of L p relative to the principal axes of I cr . The full set of fixed points, and their 
possible observable significance, will be discussed completely elsewhere; here we focus on the particularly 
simple - but far from general - situation in which L p is along one of the principal axes of I cr . In that 
circumstance, the fixed point solution is f2 cr = J~i c = fi, with J7 || L p . 

Perturbations about this fixed point are studied most easily in the frame corotating with the crust, 
where I cr is independent of time. For definiteness, let us suppose that the principal moments of inertia of 
the crust are I\ < I2 < I3 and, to parallel Shaham (1977) as closely as possible, suppose that the fixed point 
corresponds to rotation about the 3-axis at angular velocity ft = e^fl. Then the linearized equations are 

K 



n 



cr,2 









v 




h 










1^ 






h 





f^cr,2 — _ 77;(^cr.l — ^c,l) 

K 



n 



cr.l 



(n cr>2 — fi c ,2) 



^cr,3 r o (^cr.3 ^c,3) 

-<3" 

K 

^c,l + Qcr,2 — Qc,2 = "T7r(^cr,l _ ^c,l) 



K 



Cl Ci 2 + ^c,l — ^cr,l — YJj^ cr 2 ~ ^,2) 
K 

^c,3 = YJj^ cr - 3 ~ ^c,3)- 



(5) 

(6) 
(7) 
(8) 
(9) 
(10) 



In these equations, F = Q~ 1 d*F/dt, where d*F/dt is the time derivative of any vector F as seen in the 
frame rotating with the crust. It is clear that the perturbations along the 3-axis decouple from those along 
other axes, and decay exponentially with a characteristic rate 



(ii) 
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this coupling time has been estimated by Alpar and Sauls (1988), according to whom Qt/2it w 400 — 10 4 
(and the relaxation of the electron distribution function is due to the scattering off the neutron vortex 
magnetization), and by Sedrakian and Sedrakian (1995), who find that Q,t/2tt is rather sensitive to the 
mass density p and spans the range Qt/2tt <~ (10 2 — 10 8 )P(s) for p <~ (1.6 — 3) x 10 14 g cm -3 (here the 
electron scattering is off the proton vortex clusters coupled to the neutron vortex lattice). In the latter case 
the relaxation time spans a wide density range with a large gradient near the crust-core interface because 
of an exponential dependence of the size of the cluster on the proton effective mass. (See also Sedrakian et 
al 1995, Sedrakian & Cordes 1998; for the decay of precession the effective coupling rate 7 is a weighted 
average of the range found by Sedrakian and Sedrakian 1995, implying an effective coupling time closer to 
the smallest values; we adopt Qt/2tt <~ 100P(s) for numerical estimates below.) 

The remaining four equations have normal modes proportional to exp(pSlt) = cxp(p</>), where (f> = Q,t 
is pulse phase. It is straightforward to solve for the modes of a triaxial star, but the basic result can be 
derived under the assumption of axisymmetry, I 2 = h- (We have solved the corresponding triaxial problem, 
and there are no qualitatively different modes for slowly rotating neutron stars.) If we define 

h\ , L p _ J _ _K ( I ^ l \ _ h(h+I c ) 



a= {—) + m and ^ = n{T 1 + T c ) = nn 1 ( h + i c y (12) 

then 

^cr.l + C^cr,2 = _ 7 —p(^cr,l — ^c,l) (13) 

ll + i C 

lie 

f^cr,2 — C^cr,l = —7 — r"(^cr,2 — ^0,2) (14) 

^1 + i C 

7/1 

Q Cj l + Q cr .2 — ^c.2 = 7 —r(Qcr.l — ^c.l) (15) 

h + i c 

7/1 

f^c,2 + f^c,l — ^cr,l = 7 ; — p(^cr.2 _ ^0,2)- (16) 

Jl + lc 

If 7 = 0, so the crust and core are uncoupled, then there are modes with p = ±ia which correspond to 
independent precession of fl cr but at a frequency that is much larger than the conventional Euler frequency 
for reasonable values of L p /I\Q, = I p /h (where I p is the moment of inertia of pinned superfluid). The 
remaining modes with p = ±i are an artifact of working in the frame that corotates with the crust, and 
correspond to J7 C fixed in the inertial frame of reference. 

When 7^0, the modes are damped, as was discussed by Bondi and Gold (1954) in the context of 
the rotation of the Earth (without considering pinned superfluid, of course!). The characteristic equation 
is fourth order in p, but we expect the roots to come in complex conjugate pairs, so we can reduce the 
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characteristic equation to second order by introducing the complex angular velocities 

f^+ } = O cr ,i + ifl cr ,2 and f^ +) = tt cA + iO Cj2 , (17) 
which satisfy the equations 

hw-iaflw = - 7 ^ r (nw-n(+)) (is) 

1\ + lc 

- in<+> + inw = A(flW - n(+)). (19) 
Ji + j c 

Substituting oc exp(p</>) we find 

p 2 +p[7 + *(l-^)]+^(l-^^-) =0. (20) 

The normal modes of the fourth order system are the two solutions to this quadratic equation, and their 
complex conjugates. 

Although we can solve the second order characteristic equation exactly, it is more instructive to find 
approximate solutions valid for small and large crust-core coupling. For small values of 7, we rewrite the 
characteristic equation as 

{p + i){p-ia)+ 1 (p- t^V) =0; (21) 



h+l 

this form separates terms of zeroth and first order in 7 explicitly. To first order in 7, the solutions are 

_ TP + , /,/(/, +!.)] 
1 + er 



and 



^"- (ii + ZcXi + g) - (23) 



For large values of 7, we rewrite the characteristic equation as 

p-ia—!±-+ 1 - 1 \p 2 +ip(l-a)+o] =0; (24) 
n + J c 

this form is useful for expanding in powers of 7 -1 . In this case, the solutions to first order in 7" 1 are 

^ = -7-^1-7^1 ( 25 ) 



Ic + h 



and 



Pp = hTT c -WTi7){ 1 + hTTj- (26) 

In each case, pd represents damping of the angular velocity difference between crust and core, and p p is the 
precessing mode. 
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For the coupling times estimated by, for example, Alpar and Sauls (1988) or Sedrakian and Sedrakian 
(1995), the small 7 limit is the relevant one. Since I p /h ~> h/h — 1, the precession period is far smaller 
than for a rigid body, approximately h/I p spin periods. Moreover, the wobble damps away, lasting ~ 7 _1 
precession periods: 7 _1 w 400 — 10 4 according to Alpar and Sauls (1988) and a reasonable estimate for the 
effective coupling is 7 _1 ~ 100P(s) for Sedrakian and Sedrakian (1995). Even if 7 were large, the precession 
period would be short, although it would be lengthened by a factor 1 + I c /I\ relative to the small 7 limit, 
implying a cycle (I c + I\)/I p spin periods long. The precession would persist for approximately Ii/I c Qt 
precession periods in this limit. Since the crust-core coupling time must exceed the light travel time across 
the star, r > R/c s=s 0.03 ms, and the damping time for the precession must be <5000(/i/J c )P precession 
periods, where P is the rotation period in seconds. 

In neither limit is the precession either long period or persistent. From this pessimistic result, one 
concludes that free precession cannot account for the cyclical behavior seen in long time monitoring of 
some pulsars. Moreover, to explain the data, one must invoke an excitation mechanism that acts relatively 
continuously, since it must fight the tendency for neutron star wobbles to decay rapidly. The characteristic 
cycle timescales of order months to years observed for these pulsars must reflect the underlying processes 
responsible for the continuous excitations. 

2.2. Imperfect Pinning 

In demonstrating that persistent, long period precession is impossible for neutron stars with pinned 
superfluid, Shaham (1977) assumed perfect pinning. In actuality, superfluid vortex lines will not pin to 
crustal nuclei absolutely. One purpose of this paper is to see whether there are new oscillatory modes that 
emerge when pinning is assumed to be strong but not perfect. 

To study this problem, we adopt a somewhat idealized approach. In actuality, the pinning of crustal 
superfluid is a highly inhomogeneous process involving the interaction of individual vortex lines and crustal 
nuclei. This coupling is modelled by effective potentials highly localized around discrete pinning sites in 
the vortex creep picture (e.g. Anderson and Itoh 1975, Alpar et al. 1984a, Link and Epstein 1991, Link 
et al. 1993) and by scattering of particles by and Kelvon excitation of moving vortex lines not pinned 
to crustal nuclei (e.g. Epstein and Baym 1992, Jones 1991, 1992). In our calculations, we use smoothed 
hydrodynamical equations to describe the coupling between the superfluid and normal components of the 
crust macroscopically, using the formalism developed in Khalatnikov (1965, Section 16). This formulation 
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of the problem is linked most naturally to a picture in which supcrfluid vortex lines experience drag forces 
as they move through a smooth medium of normal fluid, but also may be applied directly in the vortex 
creep picture in the linear approximation (i.e. when the difference between the angular velocities of the 
superfluid and normal fluid are sufficiently small). 

Shaham's results are recovered in the limit of perfect coupling, that is, when the coefficients of mutual 
friction are infinite. We can explore whether qualitatively new modes appear when the mutual friction is 
strong, but pinning is not perfect. As we shall see, no new slowly damped, long period modes arise. 

General formulae for the mutual friction force are given in Khalatnikov (1965, Section 16). If the 
superfluid vorticity is defined to be u> = V x v s , where v s is the superfluid velocity, then the net force per 
unit volume acting on the superfluid is 

f = -u> x (V x Xv) - (3' p s u> x u - (3p s v x (u> x u) + j'p s vuj ■ u, (27) 

where p s is the superfluid mass density, v = u>/|u>|, and, for a normal fluid velocity v n , 

u = v n v s — — V x Xv- (28) 

Ps 

X = (p s k/4:7t) ln(d/£), where k is the quantum of circulation per vortex line, d is the effective intervortex 
separation, and £ is the coherence length. The mutual friction force is defined to be f + u> x (V x Xv)/p s . 
The parameters (3 and jPrime must be positive for the rate of energy dissipation resulting from f to be 
greater than zero locally. 

Qualitatively, the terms in f involving u arise from the bending of vortex lines, and shall be neglected 
here. Of the remaining contributions to f , the two proportional to (3 and (3' are perpendicular to ui whereas 
the one proportional to 7' is along us; the latter is expected to be small, and we neglect it too. With these 
simplifications, the form for f used in this paper is 

f = -P'p s u> x (v n - v s ) - (3p s v x [u; x (v n - v s )]. (29) 

For getting a qualitative feeling for the relative sizes of the phenomenological quantities (3 and (3' , we use 
a different parametrization for the strength of the mutual friction force, based on the idea of vortex drag. 
The equation for the superfluid velocity including mutual friction is 

-£+v.-Vv a = -V(A* + 0) + -, (30) 

where \x is the chemical potential and <fi the gravitational potential; taking the curl of this equation gives 

^ = Vx(v s x U + f/ ft ). (31) 
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If f = 0, then the supcrfluid vortex lines comove with the superfluid, but, in general, the vortex lines have a 
different velocity, vl 7^ v s , and 

— = V x (v L x w); (32) 
from the form for f given in equation (p9|), we can read off 

v L = v s + /3'(v n - v.) + fiv x (v n - v.) = v n + ((3' - l)(v n - v.) + f3v x (v n - v.). (33) 

Only the components of v n — v s perpendicular to v contribute to vl, as can be seen from the original 
expression for f. Clearly, vortex lines comove with the superfluid if |/3'| and (3 are both small, and comove 
with the normal fluid if \(3' — 1| and (3 are small. (If superfluid rotates faster than normal fluid, vortices 
move slowly outward relative to normal fluid for [3 small.) The motion of a vortex is found by balancing the 
Magnus force due to superfluid streaming past the line and any other forces it experiences; for our purposes, 
the latter are drag forces perpendicular to the line, so the equation of motion is 

p s KV x (v L - v s ) + F d = 0. (34) 

If the drag force per length on a vortex is 

F d = -f)(v L - v„) - T)V x (v L - v„), (35) 

then the vortex line velocity is 

vl = v + h 2 -V(p s K-?7 , )](v n -v s ) + T)Kp s V x (v n - v s ) 
(p s K — rj') 2 + rf- {p s K — T]') 2 + rf ' 



13' = 1- (38) 



from which we infer the relations 

(p s K - if ) 2 + 7] 

p s n{p s n - f]') 
(p s K - ri') 2 + 77 

These results relate the drag coefficients 77 and 77' with the parameters (3 and (3' appearing in the mutual 
friction force. 

In microscopic models for mutual friction developed so far, the coefficient 77', which determines the 
magnitude of the drag force perpendicular to the motion of a vortex line through the normal fluid, is 
negligible. If ?/ = 0, then equations (^) and ( |38| ) simplify to 

13 = (Ps^lv 2 ' m 
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= 1- , {p :f 2 . (40) 

From these relationships, we find that when r\ 3> p s n vortex lines are dragged effectively, and tend to 
follow the normal fluid closely; in that limit, (3 Rj p s n/rj and 1 — j3' Rj (p s K/r]) 2 Rj /3 2 . When 77 <C /3 s k the 
drag is weak, vortex lines tend to follow the superfluid, and /3 Rj 77/ p 8 K and /3' Rj (r\f p s n) 2 Rj /3 2 . As we 
shall see below, this means that the dissipative torque arising from mutual friction is much larger than the 
nondissipative torque when the drag is either very strong or very weak; the two torques are only comparable 
when r// p s n ~ 1. When 77' 7^ 0, the situation becomes more complicated, as equations ( [3?] ) and ( |38| ) involve 
two nondimcnsional parameters, r\j ' p s j 'k and 77' / p s n. If we assume that rf <c 77, then in the strong damping 
limit, j3 w p s K,/r] as before, and 1 — /3' Rj (p s n/ri) 2 if 77' -C p s K, but 1 — /3' Rj —rj'p s K/r] Rj — (r)'/rj)[3. In 
the weakly coupled domain, /3 Rj r\j p s n as before, but /3' Rj /3 2 only when 77'/^ <C r)/p s K, and instead 
w —r)'/p s K Rj —(77' when 77'/ 77 3> /O s k. As we shall see, these results will make the existence of 
long term oscillatory modes problematic when 77' <C 77 provided that at least part of the crust is coupled 
strongly to the crustal superfluid. For 7/ > 77, the situation turns out to be more favorable for the survival 
of oscillatory modes. In that case, (3 Rj rjp s n/ri' 2 and 1 — f3' Rj —p s n/ri' Rj — (rj' /rj)/3 in the limit of strong 
coupling, and (3 ~ rj/p s K and (3' Rj —rj'/p s K Rj —(rf/rj)f3 in the limit of weak coupling. 

The torque that results from mutual friction is 

N = J d 3 rr x f(r) = N + Np>, (41) 

where, from equation (p9[), 

Np = - J d 3 r (3 Ps r x{i/x[wx (v n - v.)]} (42) 

N P' = - J d ^P'Ps* x [w x (v n - v s )]. (43) 

We restrict ourselves to a uniformly rotating normal fluid, but the analogous restriction to uniformly 
rotating superfluid is dynamically inconsistent unless (3 and [3 1 are independent of position. Consequently, 
we imagine that the star can be divided into "shells" in which (3 and [3 1 are independent of position, and the 
superfluid rotates uniformly. In these shells, v s = fl s x r and v n = f2„ x r, with tt n and fl s independent 
of r; this also implies that u) = 2tt s . In succeeding sections, we consider stars with one and two superfluid 
shells. These examples suffice to illustrate the complex behavior that may arise in a real neutron star, 
where (3 and (3' vary continuously. 

For uniform rotation, equations ( ff2| ) and ( |43| ) become 

N^fi s xT r (n a x n„) + n s (n s - n n ) • [T fi - stt(t )] (44) 
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N/3' = —(Sis - Sl n ) X Tpi ■ Sl s , (45) 

where 

T f3 = 2/3 J d 3 rp s rr, (46) 



-/8' 



= 2/3' / d 3 rp s rr, (47) 



5 is the unit tensor, and Tr(Tp) is the trace of T^g. It is easy to show that (Sl s — Sl n ) ■ = and 
(Sl s — Sl n ) ■ < 0, so that (Sl s — Sl n ) • N < 0, that is, mutual friction torques are ultimately dissipative. 

The tensors Tp and Tp/ can be rather complicated in general. Even in uniformly rotating superfluid 
shells, the superfluid density p s (r) may be slightly anisotropic, principally as a result of rotational flattening 
perpendicular to Sl s , which is time varying and not aligned with any of the principal axes of the crust in 
general. However, we shall neglect these complications, since the magnitudes of the anisotropics in Tp and 
Tpi are expected to be small for slowly rotating neutron stars, which we focus on here. Accordingly, we 
approximate 

T p = I s [3 cS 8 (48) 

Tp, = I s p' cS S, (49) 

where I s is the moment of inertia of the superfluid, and /3 e ff and (3' eS are suitably averaged /3 and /?'; 
henceforth, we drop the subscript "eff." With these expressions for and Tpi the mutual friction torques 
simplify to 

N p = -i s [3n s {n s -n„)-{5 + n s si s ) (50) 

%=//(fi„xfl s ). (51) 

We shall devote much of the remainder of this paper to examining the consequences of torques of this form. 
Here, we neglect other torques which could be important, such as gravitational torques (both Newtonian and 
post-Newtonian) or fluid torques arising from boundary conditions, and ignore the various complications in 
Tp and Tp> alluded to above. Some of these issues will be discussed in a subsequent publication. 



3. TWO COMPONENT STAR 



Implicit in the review of Shaham (1977) presented in Section 2.1 was a treatment of the two component 
system consisting of the rigid crust and pinned crustal superfluid. This was the 7 = limit in which 
the crust and core decouple entirely. In that case, we found that the crust precesses at a frequency a 
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(see eq. under the additional assumption of axisymmetry; for 7 = 0, this mode is undamped. The 



three component model discussed in Section 2.1 also reduces to a two component system when 7 — > 00, in 
which case the core and crust are coupled perfectly, and must corotate. Although this limit is not realistic 



(see discussion of the maximum possible 7 in the penultimate paragraph of Section 2T), it also leads to 
undamped precession at a frequency al c j(l\ + I c ). Here, we examine how imperfect pinning alters these 
results, and introduces both damping and new modes. 



3.1. Free Precession Reexamined 

The coupled equations for the angular momenta of the crust and crustal superfluid are 



dt 



= -N/3 - n i3 , = i s pn s (n s - n cr ) ■ (8 + n s n s ) + i s p'[si s x n cr ) (52) 
n + Npi = -i s pn s (n s - n cr ) ■ (s + h s n s ) + i s p'(n cr x n.), (53) 



s dt 

where we have substituted fl cr for Q n in equations ( |5C| ) and (|5l|), and assumed that the angular momentum 
of the superfluid is I s ^t s , which is tantamount to assuming that the moment of inertia tensor of the 
superfluid is of the form 

I s = (I. - AI S )S + AI S ( 3 "'"'~* ) . (54) 
Notice that if ft' = 1 and (3 = 0, these equations reduce to 

d(Icr ;" cr) +J s (n C rxn s ) = o, (55) 

dt 

which is equivalent to equation ([l]) with the contribution from the core component omitted, and 

= si cr x n si (56) 

which implies that £l s is fixed in the reference frame that rotates with the superfluid. This is the limit of 
perfect pinning, and results in undamped precession at the frequency <r. 

When ft' ^ 1 and P ^ 0, there are additional modes. Let us work in the frame rotating with the crust, 
in which case (recall that d*F/dt is the time derivative of F in this frame) 

. + o cr x (i cr ■ n cr ) = i a pn s (n s - n cr ) ■ (s + n s n s ) + i s p'(n s x n cr ) (57) 

dt 

+ (1 - p'){n cr x n.) = -pn s {n s - n cr ) ■ (s + n a n a ). (58) 
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If we project equations ( p7| ) and ( pq ) along the principal axes of the crust, and linearize around the fixed 
point at which f2 cr = £l s = $1 and fi || e"3 we find 



cr,2 



h 
h-h 



h 



^cr.2 -j—Qs.2 



(f2 crj i — ri s ,i) 



—f— ^cr,l H 7— ^s,l — J—(^cr,2 ~ ^3,2) 



^cr,3 ^ (^cr,3 ^s,3) 

^3 

- (1 - /3')(" S ,2 - fW) = " Ocr.l) 

O s ,2 + (1 ~ - ncr,l) = ~/3(^,2 " ^cr,2) 

Qs,3 — — 2/3(Sl Sj 3 — Q c r,3)' 



(59) 

(60) 

(61) 
(62) 
(63) 
(64) 



As before, the evolution of the perturbations along the 3-axis decouple from those along the other axes, 
and decay exponentially; the rate of decay is 2/3(1 + I s / The remaining equations imply a fourth order 
characteristic equation if we search for modes cx exp(p</>). 



3.1.1. Axisymmetric Crust 



When the crust is axisymmetric, I\ = I2, and equations fl59|), and (60) simplify to 

'h-h\ , h/3' 



h J h 
h-h\ , hp' 



h ) h 

these couple to equations (|62j) and (p3[), which are unchanged. 



^cr.2 -j—Qs.2 — j—(&cr,l ~ ^s,l) 

1\ ^1 

^cr,l H J— fis,l = 7— (^cr,2 — ^8,2); 

Jl ^1 



(65) 
(66) 



As we found in Section 2^ , the fourth order characteristic equation may be reduced to second order in 
this case. Define 

'h-h\ , hp' 



a 



hP' 
h 



and let 



then we get the two coupled equations 



n { s +) = n sl + m 



s,2i 



(iw - io>nw + i^n« = -^(nw - nW) 



(67) 



(68) 



(69) 
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+ i( i _ /3')(oi+) - 4+)) = -/?(n<+> - nW). 

It turns out to be convenient to use 



instead of r^i -1- ' ; doing so yields the coupled equations 

SI™ - iefi<+> + ^(^' _ ^)A(+) - 



and 



A(+) 



Assuming that (Ocr , A' + )) cx exp(p0) we find the relation 



A« = 



(+) 



p + tfl-jS'J+jS' 



and the characteristic equation 



P + P 



P' -a')+p 1 + 



+ e(l - /3' - i/3) = 0. 



(70) 



(71) 



(72) 



(73) 



(74) 



(75) 



Equation (|74|) is useful for finding the eigenvectors once equation (|7^) is solved; these are needed to 
determine the response of the two spin components to external torques. 

Although we can solve equation (|75| ) exactly, it is instructive to consider the two limiting cases of weak 
and strong vortex drag separately. When vortex drag is weak, (3 and (3' are small in magnitude, so we 
rewrite equation J7^ ) as 



X, 



ip[ 1 + -^ ) +e 



0. 



The solutions to this equation to first order in the small quantities (3 and f3' are 

(ip-P)(l + e + I a /Ii) 



and (the Shaham mode) 



Pd = -i + 



P P 



1 + 6 



e(ip' - (3)I S 



(76) 



(77) 



(78) 



Both of these solutions damp slowly, at rates proportional to f3 > 0. The second mode reduces to the 
conventional Euler precession when /3 = j3' = 0. The first mode arises because the superfluid angular 
velocity would remain fixed in the inertial frame if (3 and (3' were zero, but wanders slowly when the 
coupling is small but nonzero. (We discuss this point more fully in the context of three component models 
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see Section [4.1.1| below, discussion following eq. [14S ] .) For strong vortex drag, 1 — (3' and (3 are small, and 
we rewrite equation (fTq) in the form 



Pip^ia 1 ) + (1 -f3')(ip + e) + (3 



pi 1 



h 



0. 



The solutions to this equation to first order in the small quantities (3 and 1 — (3' are 

e[(3 + i(l-(3')} 



Pd 



and 



Pp 



ho' 



(79) 



(80) 



(81) 



The first of these solutions represents a slowly damped mode with an oscillatory part that is negligible 
when T]' <C r], implying 1 — j3' <§; (3 in this limit. (Recall discussion following eqs. [M and |38 in Section 



2.2.) The second mode corresponds to precession at a' with slow damping: For I s /I\ ^> e, equation (67) 
implies that a' ~ /s/^i for 1 — /3' <C 1, so the damping rate is approximately (3(1 + I s /I\)\ in the unlikely 
event that e 3> I s /Ii, then a' ~ e and the damping rate is approximately (I s / 'I\e)(3(l + e). 



3.1.2. Non-Axisymmetric Crust 

Since the neutron star crust may not be axisymmetric, it is worth checking that there are no surprises 
when 1\ =/= h. Define 

then equations (|59|) and (p0|) become 



^cr,l + Ci^cr,2 f— ^s,2 — 1 - (^cr,l ~ ^s.l) (84) 

&cr,2 — o' 2 flcr,l H r~ ^s.l = T~ (^cr.2 — ^s.2), (85) 

J2 H 

which must be solved along with equations (|6^ ) and (|63|). When we look for solutions cx exp(p0) we find 
the fourth order characteristic equation 

V{ [/? + (i - «*] (1 + £) (1 + 1) - (i - fff^ 2 - (i - /?) (| + 1) + ✓,✓,} 

+P/3(2e ie2 + ^2 + ^2.) + Wl \S 2 + (1 - ,3') 2 ]. (86) 
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It is not possible to factorize the characteristic equation into the product of two second order equations 
because there is no guarantee that all of the roots are simply complex conjugate pairs. 

In the limit of weak coupling, expanding equation (R6|) up to first order in and 0' yields 







(p 2 + l)(;p 2 + eie 2 ) + 



p A \ 2+^ + ^ I +p[ 2eie 2 



-0'W 



h h 
7 s (l-e 2 ) , 1.(1 -ei) 



h h 



2eie 2 



h h 

The approximate solutions of this form for the characteristic equation are 

[2(1 - eie 2 ) + (7./Ji)(l - e 2 ) + (7 S /7 2 )(1 - ei)](i^ - /?) 



(87) 



and 



Pp 



V e l e 2 



2(l-eie 2 ) 
i/3'V^[(7,/7i)(l - e 2 ) + (7 S // 2 )(1 - e x )] 



2(l-eie 2 ) 
/3[(7 s /7!)e 2 (l - £l ) + (7 s /7 2 )ex(l - e 2 )] 



(89) 



2(l-eie 2 ) 

and their complex conjugates. To get limiting results in the strongly coupled domain, we rewrite equation 
(Em in the slightly modified form 







pV + <^ 2 )+p73(2 + -^ + 



I, 



V{ r + (i - 0'?\ (i + 1) (i + 1) - (i - ^ - (i - /?') (I + f ) 

-HPjflf 26 l£2 + ^ + ^ + ei£2 [/3 2 + (1 - /3') 2 ]. 



(90) 



From equation (|9C|), it is evident that the modes are near p 2 = and p 2 = —o' x a' 2 . To get the first order 
approximation to the modes with p 2 = to zeroth order in and 1 — 0', we identify all terms in equation 
( |90| ) that are potentially second order in small quantities; this leads to the quadratic equation 

7 s e 2 ( I s e 1 



Q=p 2 cr' 1 <7 2 +p0[ 2eie 2 



+ e l e 2 [0 2 + {\-0')% 



which has the pair of roots 



±_ 0(^2 + £ 2CTi) 



2a' lCT2 



± 



2 (e 1 a' 2 -e 2 a[) 2 6^(1-0') 



IY2 



{2a[a' 2 ) 2 



1/2 



(91) 



(92) 



In the axisymmetric limit, these two roots reduce to the damped mode found in Section 3.1.1 and its 
complex conjugate, but although p^ both imply damping in general, they could be purely real and different 
in magnitude, especially since we expect 1 — 0' to be much smaller than in the strongly coupled regime. 
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(This is why we could not factor eq. |86| into two quadratic equations.) It is also straightforward to expand 



equation ( p0[ ) around the approximate root p ~ i \J to find 



Pp 



to first order in the small quantities (3 and 1 — /?' 



(93) 



From this brief foray into the modes of a triaxial star, we conclude that deviations from axisymmetry 
do not alter the behavior of the precession qualitatively in either limit. The character of the damped modes 
may be different in the strong coupling limit for nonaxisymmetric stars, but if so, they become purely 
damped, with no oscillation at all. Consequently, from here on we specialize to axisymmetric crusts, since 
we do not expect to miss any important oscillatory modes. 



3.2. Response to External Torques 

As we have seen, the modes of free precession for this two component model are rapidly oscillating 
and/or damped. Here, we consider the response of the system to external torques. Our analysis will reveal 
that the limit of perfect coupling between the normal fluid and superfluid must be taken with care when 
external torques act. 

If we suppose that the crust is subject to an arbitrary time-dependent torque, N cr (</>), then equations 
( f72| ) and (|7^) are changed to 

n<+> - iert+ ] + ^{0 - /3)A<+> = N£\<f) (94) 
h 



and 

A(+) + /i[l-/3'(l + | 



+ / ?(i + |)}a(+) +ie fi(+) = -A>(+)(0), 



(95) 



where Ncr — I\ 1 (N cr j + iN C r,2)- Apart from decaying transients, the solution to these equations is 

ft(+)= Y, A * I d<f>'N^ (</>') expfr^ -</>')} (96) 

_ J J — OO 



a—p,d 



a—p,d 

where the coefficients arc 



a(+) - g Pa + n%) + L ^ eMp ^ - * % (97) 



[p p + i(l -/?)+/?] [ Pd + i(l-/3 , )+/?] 

A p = r A d = . (98) 

(Pp ~ Pd) (Pd ~ Pp) 
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(Equation [J74j is useful in obtaining the A a .) 

In the limit of perfect coupling between the crust and superfluid, we take i(l — 0') + [3 — > with 
Pd/[i(l — 0') + 0\ = —e/cr' constant; taking this limit of equations ( p6[ ) and ( p7| ) naively yields 

fi« = d(f>'N^ exp[p p (^ - 0')] (99) 

J — DO 

with A(+) -> 

Actually, the perfect coupling limit is a bit more subtle than the manipulations leading to equation 
(p9|). To see why, consider the response to a time-independent torque on the crust. Then equation ( |96| ) may 
be integrated readily and we find 

7-/V ( + ) 

= (ioo) 

integrating equation ( |97| ) yields A^ + ) = 0. These results ought to hold for any i(l — [3') + (3. However, 
equation (f)9|), which purports to describe the limit of perfect coupling, i(l — (3 1 ) +(3 = 0, yields 

/v (+) 

fiW = -— - (101) 

J3p cr' 

(assuming a small, negative real part to p p ). Which of these results is correct? 

To resolve the conundrum, consider a torque that turns on (or can be regarded as constant since) some 
time in the past, <j>Q. Then equations (|96|) and (|)7|) yield the response 

i ( i(l-/3')+/3\expbj>(4>-^o)] 



e V P P / Pp-Pd 



V Pd / Pp ~ Pd J 



A (+) = f " exp[p p (</> - <fo)] + exp[p d (0 - <fo)] 1 ^ (+) 



(103) 

Pp-Pd J 

We suppose that the damping constant associated with the precessing mode, p p , is large enough that 
exp[p p ((f> — 0o)] ~~ * 0; then what we find depends on Pd(4> — 4>o)- (Recall that the real parts of p p and pd 
are negative.) If |Re[p,j(0 — 0o)]| 3> lj so that any transient response has plenty of time to damp away 



between <j>$ and 0, then we recover equation (10C), and also find that A^ + ) — > 0. On the other hand, if 
|Re[pd(</> — 0o)]| "C 1, then we recover equation (101) and A' + ) — > — 57cr in the limit of perfect coupling, 
i(l — (3') + (3 — > 0, to zeroth order in pd(0 — <fio)- To first order, we also find a term that grows linearly 
with Pd{4> ~ <Po )'■> over a sufficiently long timespan, this growth would change the tilt from equation (|100|) to 



equation (101). 

The resolution of the apparent paradox is that there is none: equations d96| ) and (|97j) are always the 
correct ones to use. What one gets in the limit of strong vortex drag depends on how the timescale on 
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which the external torque changes compares with the timescales inherent in the coupling of supcrfluid to 



the normal crust. Equations (100) and (101) both have domains of validity; equation (101) is a lower bound 
to the steady-state tilt of the rotational angular velocity away from the 3-axis in the strong coupling limit. 
As long as the damping timescale associated with pd is short compared to any timescale associated with 
changes in the external torque, however, equation (|l00| ) gives the right response. In practical terms, pulsar 
spindown provides a nearly constant torque on the crust which can give rise to Ncr . Thus, if pd implies 



decay on timescales smaller than the pulsar spindown time, then equation (100) describes the response of 
the crust, even in the strong pinning limit. 

We emphasize that this result could not be found from a consideration of normal modes alone; arriving 
at it requires examining the response of the star to a torque. There is therefore a subtle aspect to the 
limiting case considered by Shaham (1977): While the modal frequencies he derived are correct, and his 
conclusions about free precession warranted, blithely using equation (|99|), which would follow from the 



assumption of perfect pinning, instead of equations (|96|) and (97) is wrong even in the limit of strong 
pinning. 

The response of the star to torques along the 3-axis is found by solving the equations 



cr,3 



and 

where N cr ^ 3 = N cn3 // 3 ; the result is 

A 3 = - 

and 



iv cr>3 (</>) 



A 3 + 2/?( 1 + j )A 3 = -N crt3 {0) 



t>'N cri3 ((j)') exp 



-2/3 1 



h 



'NcrAt') + T-hr / WNcrM') eX P 
J.s + J 3 . /_™ 



-2/?(l + ^ 



In particular, a time-independent torque results in a steady angular velocity difference 



2/3(1 +I a /hY 



while the crustal angular velocity changes linearly: 



N, 



cr, 3 -i 3 



(104) 



(105) 



(106) 



(107) 



(108) 



(109) 
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The response to an impulsive torque iV cr) 3 = No, 3 <5(0 — 4>o) is 

O = No, 3 {/ 3 + Is cxp[-2/?(l + I s /I 3 )(4> - <j> )}} 
cr ' 3 h+h 1 ' 

A 3 = -N ,3 exp[- 2/3(1 + J,/I 3 )fo - <&,)]• (HI) 

Spindown torques, which only change on very long timcscalcs for all observed pulsars, provide a 
physical realization of a "time-independent" torque. Let us take the vacuum magnetic dipole torque, which 
is proportional to — /j, x (u> x fi) (e.g. Davis and Goldstein 1970, Goldreich 1970, Michel 1991); assuming a 
magnetic dipole moment fi that is fixed in the frame of the crust, with components 

fi — cos a§3 + sin ae± , (112) 

the torque may be written (to a good first approximation) as 

Nd = Na sina(ei cos a — e 3 sin a) = N s d(ei cot a — e 3 ), (H3) 

implying Ncr — N sc i cot a/I\ and ^V cr , 3 = ~N sc i/l3- The steady state response to these torques is 



OW = and A<+) = 



(114) 



and 



2fj(I s + i 3 ) 



with f2 cr . ; 3 = —N s d/l3. Two aspects of these results are especially noteworthy. First, as is well-known, 
the superfluid rotates faster than the crust as a consequence of the spindown torque. Second, but not so 
widely appreciated, the steady state "tilt" in the angular velocity of the crust is surprisingly large, as it is 
proportional to (T 1 . (When e = 0, eq. |7jJ implies that pd — 0, and flit^ grows linearly with time according 
toeq. @.) 

3.3. Free Precession Encore 

3.3.1. Different Crust and Superfluid Angular Velocities in Unperturbed State 

The fact that time-independent spindown of the crust implies a difference between il cr , 3 and D Sj 3 in 
steady state suggests that we explore free precession once again, but with a different unperturbed state 



than we used in Section 3.1. There, we assumed that the undisturbed star rotates with Q s = f! cr = ^§3. 
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Let us consider instead what happens when Q c = fles and 1~2 S = £flea in the unperturbed state. We shall 
not assume that |£ — 1| must be small, although we expect this to be true. 

The equations governing the precession of an axisymmetric star given this new unperturbed state are 
0« - ia'^ip + = -^(n« - (116) 

n<+> + i(i - /?')(^ +) - ^+ } ) = - (117) 

where 



(118) 



These equations yield the characteristic equation 



P 2 +P 



-a' i ) + p(\ + I f\ +£ (l-[3'-i(3)+ i^(l-t) = 0; (119) 



in the weakly coupled limit, the solutions are 



and in the strong coupling limit, 



, (ig -/?)(i + e + er s /ii) 

1 + e 

P (mi 



jeCW0+^[e + / S (g-1)/J 1 ] 

Pd = (122) 

I s K(l-/3')+/3(l + ^)] 



These equations are hardly different than their £ = 1 counterparts except for two noteworthy differences. 



Vv = io\ —, (123) 

il<7. 



First, in the weakly coupled limit, the real part of p p can become positive for £ — 1 > e, implying linear 
instability. However, this is not worrisome since the growth time of the mode is of order the spindown time 
of the star if £ — 1 = N s d/2f3(I s + Second, in the strongly coupled limit, the real part of pd is enhanced 
for £ — 1 > 0, implying faster damping. Otherwise, the effect of £ ^ 1 is merely to renormalize the various 
coefficients appearing in the solutions to the characteristic equation without introducing any qualitatively 
new behavior. Consequently, we shall not consider the effects of differential rotation further in this paper. 



3.3.2. Interaction of Two Regions with Simple Modes 



In Sections |3.1.1| and 3.1.2, we explored the characteristic modes of precession for axisymmetric and 
nonaxisymmetric crusts in the limits of weak and strong coupling to the crustal superfluid. In general, we 
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can imagine that there are regions of the crust to which the superfluid couples with different strengths. 
These regions are presumably linked to one another by elastic forces that seek to enforce corotation. 

To get a crude idea of what might transpire in such a situation, let us imagine dividing the crust into 
two components, a and b. To get a schematic feeling for the normal modes of the coupled system, suppose 
the angular velocities can be described by the equations 

n<+> - Pa nM = -^^(n<+) - o[ + ») (124) 
^"^ +, = "^ +) " ni+)) ' (125) 

where Sl a = fi ,i + i^b,i an d ^j, = + ^6,2- This system of equations has the characteristic 



equation 



= P -P{Pa+Pb-V+PaPb T— -j • ( 126 ) 

l a + h 



When 7 is small, we rewrite the above equation as 



() = (/; [p ; 2'' ) w 



to first order in 7, the roots are 

PI = Pa" (128) 

P2= Pb - T ^ T . (129) 

The effect of weak coupling is additional damping of the modes of the individual components. When 7 is 
large, we rewrite the characteristic equation as 

= p Pja+ Pbh +1 -\p- Pa )(p Pb ). (130) 

la + lb 

In this limit, the two roots are 

. Pah+Pbh no1 . 

Pi = -7 + j.j (131) 

la + lb 

Pala + Pbh (Pb - Pa)(Ia - h) H on s 

P2 = -laTTb Wa + h) ■ ( ] 

In this case, p\ represents almost pure damping at a rate close to 7, whereas P2 represents damped precession 
at a rate that is approximately the sum of the precession frequencies of the individual components weighted 
by their moment of inertia fractions. 

Consequently, if the characteristic coupling timescales among different components of the crust are 
short, we find a mean precession frequency weighted toward regions that comprise the bulk of the crustal 
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moment of inertia. Only if the coupling times are long will the precession frequencies of individual crustal 
components be apparent. This indicates that even if there are small regions of unpinned superfluid in the 
crust, precession at the Euler rate appropriate to those zones may only be seen if they do not couple to 
the rest of the crust efficiently, and that even if this is so, the precession will damp out eventually as a 
consequence of the dissipative interaction. 

From time to time, it has also been suggested that Tkachenko modes of the core superfluid could 
manifest themselves in the observed rotation of a pulsar, which is presumably the angular velocity of 
its crust. If the crustal superfluid is strongly pinned, we can regard one component, a, as the crust, 
with p a = a = e + I s /I±, and the other as the core, with p the oscillation frequency of the core in the 
zero-coupling limit. Equation ( |129| ) shows that if the interaction between crust and core is weak, there 
is indeed a mode of oscillation close to pb- For small but nonzero 7, the oscillations decay at a rate 
jl a /(h + h), which implies a damping time ~ 7~ 1 / c /7i rotation periods, if we assume I a = I\ Ij, = I c ; 
this would be of order (400 — 10 4 )/ c //i for the damping times estimated by Alpar and Sauls (1988), or 
approximately days to months for I c /Ii ~ 10 2 , comparable to the estimated periods of oscillation if they 
can occur. Moreover, for small values of 7, the effect of the core oscillations on the crustal angular velocity 



is small: equation ( |124| ) implies that 



nl+) = 



lh 



n[\ + \ (133) 



Ah + Ia){Pb - Pa) + lh 

which means that the amplitude of the oscillations in the angular velocity of the crust is ~ jh/Pa(h + la) 
times the amplitude of the oscillations of the angular velocity of the core, assuming that \p a \ ^> \pb\- Thus, 
even if the coupling is so weak that the damping time td associated with core oscillations is very long, their 
observable manifestation in the angular velocity of the crust is suppressed by a factor ~ (Cltd) . 



4. THREE COMPONENT STAR 

The calculations presented in Section || are valid if the crust and crustal superfluid are either completely 
decoupled from the core of the star or coupled to it perfectly. In the limit of perfect pinning, Shaham (1977) 
found that precession was damped as a consequence of imperfect coupling to the core. Here, we examine a 
model consisting of three components, two of which are superfluid components that couple directly to the 
rigid crust. We have in mind two possible applications, one in which the three components are rigid crust, 
crustal superfluid and core (super)fluid, and another in which the three components are rigid crust and two 
different regions of crustal superfluid with different frictional couplings to the rigid crust. 
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4.1. Free Precession 



4.1.1. Crust, Crustal Superfluid and Core (Super)fluid 

We assume that the core couples directly only to the rigid component of the crust, via a torque of the 
form 

n cc = ~i c (n c (n c - n cr ) + i c C(n cr x n c ), (134) 

for small differences between the angular velocities of the crust and core, which are both nearly fle^. This 
form of the torque is analogous to equations ( |50|) and (|5l|), except that the dissipative torque has been 
assumed to be isotropic (by contrast to N/3). We have also included a non-dissipative torque in N cc , unlike 
Shaham (1977; see Section 2.1); this contribution can be ignored by setting £' = 0. Below, we shall assume 
that C ~ C'j a t least for keeping track of small quantities. 



If we define 

then the coupled equations for the three component star may be written in the form 



A (+) 4 < / 

A1+) + < ! 



i-C 1 



- 0(1 + 1) Jaw + ienw _ | (lC ' _ C)A (+) = 

I)] + C(l + |) }AW + idlW - T p(3' - /3)A« = 0. 



Assuming modes proportional to exp(p</>), we find that 



A<+' = - 



pQ 



and A< +) = -- 



(+) 



p + i(l - /3') + /? "" >l " c P + i(l-C0 + C 
The third order characteristic equation for this system may be written in the form 



= {p + i)^p 2 +p 0' -a')+0\l + j^J 



+ e(l -0' - 10) 



-ac-c) 



p 2 1 + -M +p z 



il 



+ /3 1 + 



h 



+e(l-0'-i0) 



(135) 

(136) 

(137) 
(138) 

(139) 



(140) 



where a' is defined in equation (67) 
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When ( and (' are small, two of the solutions of equation (14C) are close to the solutions pd and p p of 
the second order equation ([75]); to first order in £ and the corrections to pd and p p are 

(I c /h)(iC - 0\Pd(icr' - Pls/h) - e(l i/J)] 



Spd 



Sp p 



{i + Pd)(pd - p P ) 
(I c /h)(iC - C)M^' - Pl./h) - e(l - (3' - i(3)} 



(i + Pp)(P P ~ Pd) 

respectively. The third solution is a new mode near p = — i; to first order in small quantities it is 

- - C)[l + e + (I s + I c )/h] 



Pd 



-1 + 



(i+Pd){i+P P ) 

When (3 and ft are small, we use equations ([77]) and (|7|) in equations ( |l4l|) and (142) to find 

UK' -Q{e + I s /h) 



5pd 



the new mode is 



Pd 



-i + 



h{l + e) 
I c e(iC - C) 

6pp = T(TTT ; 



(iC'-Q(l + e + I a /I 1 +I c /I 1 ) 



(141) 
(142) 

(143) 

(144) 
(145) 

(146) 



l + e + I 5 /Ji 

When j3 and 1 — f3' are small, upon using equations (^o| ) and (|8l| ) in equations (141) and (142) we find that 
Spd is higher order in (3 and 1 — 0', and hence very small, while 



Sp p 



Pd = -I 



h{l + a') 
(*C- C)(l + <r' + I c /h) 



1 + <T> 



(147) 



(148) 



these two results are equivalent to equations ( P3| ) and ([22j) of Section 2J if we substitute a for a' and —7 
for (i£' — £)(! + Ic/Ii)- The qualitative conclusion reached by Shaham (1977) for weak crust-core coupling 



is duplicated here: according to equations ( |145| ) and (147), precession damps out in ~ h/I c C precession 
periods irrespective of the effectiveness of vortex drag. We note, though, that the mode corresponding to 
p' d implies angular velocities that are nearly but not precisely fixed in the inertial frame of the observer; if 
(' > these could complete at least one period of oscillation before decaying away (although we do not 
expect this to be the case generally). Such modes are also found in studies of the rotation of the Earth, 
where they may arise from departures from rigid rotation of the fluid core confined by the overlying crust; 
for the Earth, the result is a retrograde motion of the pole (see Lambeck 1980, Section 3.3 for a physical 
and historical review). Since the mode arises because the core angular velocity remains fixed in the inertial 
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frame when £ and f are zero identically, we expect that for small but finite crust-core coupling, the mode 



corresponds principally to oscillations of the angular velocity of the core. From equation (139) we find 



( + ) 



(C+ + 1) « ( c+ic)(^^) 



(149) 



for this mode, which decreases linearly with the frequency of oscillation, but may be substantial nevertheless 
if/ s + 7 c >J 1 . 



(150) 



For large values of |if — f|, one of the roots of equation (|140| ) is 



Pd = (K'-C)(i + ^ 



i 1 



ic + /l , 

which is equivalent to equation (^5|) if we substitute a for a' and —7 for (if — £)(1 + I c /I\); this root does 
not depend on the strength of the vortex drag to lowest order in |i£' — CI -1 - When vortex drag is weak, the 
other two roots of equation ( |140| ) arc 

{0 -0)(1 + e + I a /h +h/h) 



P P = 



1 + e + 4//! 
I s e{i(3'-(3) I c e(i? + + e + IJ h) 



(151) 



(152) 



1 + I c /h h(l + e + I c /I x ) h[(C') 2 + C 2 ](l + h/hf ' 
(Corrections to eq. [151] proportional to |if — are higher order in (3 and f3', and have been dropped.) 



In the limit of strong vortex drag, 



Pd 



[i(l-?) + 0\ 



1 + e'- 



er' 

Ur. 



(153) 



I c cj'{iC + Q(1 + e + Ic/h) 

j 1 [(C') 2 + C 2 ](i + /c//i) 3 ' 



(154) 



(Corrections to eq. [ 1 53 1 proportional to |if — CI 1 are higher order in and 1 — /3', and have been 



dropped.) The correction for strong but imperfect crust-core coupling in equation ( 154 ) is equivalent to 
equation ( p6| ) if —7 is substituted for (if — C)(l + Ic/h)' 

Interaction with the core of the star enhances the damping of the precessing modes in all of the limiting 
cases explored above. For very small or very large coupling between the crust and crustal superfluid, the 
crust-core interactions are the principal cause of decay. In the strongly pinned regime, the characteristic 
timescales for decay are just what Shaham (1977) estimated. Imperfect pinning allows a new eigenvalue pd, 
but the associated mode damps out quickly, and so cannot be the explanation for observations of persistent 
cyclical variations in pulsar spin rates. 
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4-1-2. Crust Coupled to Two Different Regions of Crustal Superfluid 



The equations governing the spin dynamics of this system are the same as equations (136), (137) and 



(138) if we identify Jl c with the angular velocity of the second crustal superfluid component, and C arL d 
£' with the coefficients coupling this component to the rigid crust. Then it is clear that the characteristic 



equation for the normal modes of this system is still equation (140), which we rewrite in the form 







p 3 - ip 2 ( e + 



ip 2 [l + 



+ (1 - - i/3) 
+(i-C'-<) 
+(l-/3'-i/?)(i-C'-0 



■pi e 



k 

' h 

^,l + g, +P .e + | 



ie — p\ 1 



h 



(155) 



which exhibits symmetry under interchange of superfluid components explicitly. 



This form of the characteristic equation is especially useful when both superfluid components couple 
strongly to the rigid crust. In that limit, one of the roots is 



p p = 10 

to first order in small quantities, where 



(156) 



(157) 



the appearance of this root suggests a simple generalization to a multicomponent superfluid, with separate 
moments of inertia 7 SJ and coupling coefficients /3j and : 



P P 



where 



E I 8, j 
h ' 



(158) 



(159) 



The other two roots are first order small to leading order; they are approximately equal to the two roots of 
the quadratic equation 

Ir 







p 2 o~" + ip 



(l-/3'-H3) e + 



e (l-/3'- i /3)(l-C'-<)- 



(160) 
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Whcn the crust is only slightly nonspherical, soe< 1, the two roots of this equation are approximately 

- is') + moth) + [*(i - co + m/h) 



p+ 



[i(i -/?') + -C) + C] 



(161) 



(162) 



' ' [i(l - (3 1 ) + /3](I c /h) + - C) + (}{I s /h) ' 
Notice that when one of the superfluid components is coupled to the rigid crust more strongly than the 

other, p + is dominated by the tighter-coupled component, but the slow mode p_ is dominated by the 

weaker-coupled one, and reduces to equation (p0[), with a' w I s /I\ if I s \i(l — CO + CI ^ / c |i(l — + /3| and 

a' » Je/Jj if / c |i(l - /?') + 0| » J,|i(l - CO + CI- 

The limit in which one component couples strongly to the crust while the other couples only weakly 
is also of interest potentially, particularly in the aftermath of a pulsar glitch, in which some parts of the 
crustal superfluid may decouple rapidly and recouple only slowly if at all (e.g. Sedrakian 1995). The modes 



for that situation are given by equation (fell) with the correction given in equation (147), equation (pOJ) and 



equation (148). If C' > it is possible that the mode corresponding to p' d given by equation ( |148|) yields 
slowly damped oscillations in the angular velocity, as was discussed in Section [4.1.l| . 



4.2. Response to External Torques 

This section is analogous to Section ft. 2] , except that we need to consider three distinct external torques, 
acting on the crust, crustal superfluid and core, respectively, and the linear response of the three different 
components to each. 



4-2.1. Response to External Torques on the Crust 

When the crust is subject to an external torque N cr (</>), 

hip - iefiW + k{i(i> - /3)A(+) + h-W - C)A« = N£\4>) 

A<+) + j*[l - + |) ] + |) }a« + iefiW - |(iC' - C)Ai +) = -JV£f>(0) 



Ai+) 



i-C i 



h 



+ C(l + |) }Ai +) + idl™ - - l3)A (+) = -N^(<f>). 



It is straightforward to show that these equations have the particular solution 



j if J — oo 



(163) 
(164) 
(165) 

(166) 
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A(+) = ~ E „ ..n^ufl f d<t>'N£H4>')exp\p a (<t> -</>')] (167) 



A ^ +) = - E r^77vT7 /" exp[p«(4> - <t>% (168) 
where the coefficients are 

Ap = lp P + i(l-l3') + f3}lP P + ^ -(') + (} (16g) 

A d = ~— ' I -v; ^ ' ^ ,170) 







bp 


-Pd)(p P -P' d ) 


[Pd - 


H(l 


-/?' 


) + /3}\p d + i(l-C) + C] 






(Pd 


- P'd)(Pd ~ Pp) 


[PdH 




"/?') 





A ' d - - {P' d -Pd)(p' d -Pp) • (171) 

Qualitatively, these are similar to what we found for the response of a two component star, except that the 
response timescales are shorter as a consequence of crust-core coupling, which implies (perhaps significantly) 
enhanced decay of the modes. Notice that as the coupling between the rigid crust and crustal superfluid 
becomes perfect, A^ + ) — * — Ocr , an d the response simplifies to 



n&° = (P P - Pd)- 1 / d^N^Wmpj, + i(l - C') + C] exp[p p (0 - c)>')] 

-[Pa + i(l ~ + C] exp[p' d (0 - 0')]} (172) 
A< +) = (p p -p^ 1 f dcf>'NM&'){-p p exp\p p &- <j>')] +p' d eMP d ^ (173) 

J — oo 

these results are identical with equations (^6|) and (|97]) for the response of a two component star, with the 
substitution of Ac for A^ + \ p' d for p d , and i(l — (')+( for i(l — /?') + /?. However, their usefulness is 
restricted to iVcr (</>) that vary on timescales short compared with the characteristic damping time implied 
by p d , the slowest decaying mode, as was discussed in Section |3.2j . 

The response to a time independent torque on the crust is simply flit ' = iNct^ /e, just as in the two 
component case, with A^ + ) = Ac + ' = 0. These are the approximate responses obtained for a constant 
torque originating at a finite time 4>o m the past provided that p d {4> — 4>o) 3> 1. 



.£.,2. Response to Torques on the Core (Super) Fluid 



If the core is subject to a torque N c (</>) then 

fi£> - ieQ« + f (i(3' - (3)A i+) + r(iC - C)Ai +) = o 



A(+) 



1-/3' 1 + 



+ ^l + g| A (+)+ienW-|(^-0 



Al+) = 



(174) 
(175) 
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+ {* [l - C'(l + I)] + C(l + |) }Al + ' + ieilW - |(i/3' _ /3)A (+) = (176) 
where iVc = 7 ( T 1 (N c .i + iN C) 2). Apart from decaying transients, the solution to these equations is 



A c +) ^ - E , -n^n, , T ^'^ +) (0')expK(0-0')], (179) 

p,d,d' Pa +Hi U+W-oo 

where the coefficients are 

o = b P + »(i - CO + CM + j(i - CO + CM + »(i - CO + C][p p + »(i - + fl nsm 

[i(l - CO + C][*(l -13') + 13- *(1 - CO - C](ft. - Pd)(ft> - Pd) 1 j 

o _ b P + j(i - CO + CM + j(i - CO + CM + Hi - CO + CM + *(i - 10 + , , 

[*(i - CO + CMi -13') + 13- »(i - CO - C](p«i - ftOCw - pi) 1 1 

[ Pp + z(i - CO + CM + *(i - CO + CM + »(i - CO + CM + »(i - P) + / 1S21 

d [*(i - co + ci[i(i - jso + j9 - <(i - co - d(pi - ft.)(pi - j»d) ' 

When the crust and crustal superfluid are coupled to one another perfectly, pd/[i(l — (3') + (3} — > —e/a', 
and ^ -> 0, so A<+) -> -fliP and 

-exp[p^-^)]} (183) 

A ( +) = ^ + y - o + - n 



^b P + »(l-CQ + C]exp^(^-^)] 
[i(l-C) + {](p P -p' d ) 



(184) 



As was discussed in Section 3.2, this simplified version of the response may only be used for torques that 



vary rapidly compared to the damping timescale implied by pd- 



5. CONCLUSIONS 

Shaham (1977) demonstrated that when superfluid pins perfectly to crustal nuclei, the precession 
period of a neutron star is shortened immensely, and, moreover, the precession damps quickly as a result 
of weak coupling to the stellar core. One of the principal goals of this paper has been to examine whether 
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there are additional modes with long periods and long damping timescales when the assumption of perfect 
coupling between crustal nuclei and superfluid is relaxed. In fact, when the coupling is strong but imperfect, 
there are new modes that have very long characteristic timescales. One new mode is given by equation J80|) , 

for an axisymmetric, two component star, where e is the oblateness of the star and a' = e + I s /Ii, where 
I s is the moment of inertia of the crustal superfluid and I\ one of the principal moments of inertia of the 



crust. (Coupling of the crust to the stellar core hardly alters this result; see discussion in Section 4.1.) Since 
(3 <C 1, in the limit of strong vortex drag, this mode is extremely long-lived; moreover, since 1 — /?' <C 1 in 
this limit, the mode undergoes oscillations that are also extremely slow. The problem is that we expect that 
|1 — (3'\ ~ 1 <C /3 in the strong coupling domain as long as the vortex drag coefficient rj ', which governs 
the strength of the drag force perpendicular to the direction of motion of the vortex through the normal 
fluid, is small compared with rj, the analogous coefficient for the strength of the drag force antiparallel to 
the direction of motion (e.g. eqs. |37[] and |38| and ensuing discussion). Thus, this mode is not actually 
oscillatory at all, for it damps before it can complete a single cycle. In fact, for a nonaxisymmetric star, pd 
splits into two modes, with (see eq. WM) 



± _ P{^2 + £ 2gj) , 

both of which may be purely real and decaying. 



Another new mode arises in three component models when, for example, the crustal superfluid is 
coupled strongly to the rigid crust in some regions and weakly in others, or else the crustal superfluid 
is strongly coupled to the rigid crust but the superfluid core is coupled to it only weakly. Under such 
circumstances, one solution to the three component characteristic equation is equation ( |148| ) 

p , = _ i+ «--c)(iw + W) m 

where £' and £ are the coupling parameters between the rigid crust and the component that is barely tied 



to it. As was discussed in Section 4.1.1, this mode can lead to a slow wandering of the pole of the neutron 
star as seen in the inertial reference frame. However, the excitation amplitude is relatively small for the 
crustal angular velocity in this mode (e.g. eq. ]149[); moreover, in the weak coupling domain, we expect 



£' <C C if v' ^ V ( e -g- Section 2.2), so the mode decays before completing one oscillation. 

Thus, it appears likely that although there are new, possibly long-lived modes for a neutron star with 
strong but imperfect coupling between superfluid and rigid crust, these modes are not principally oscillatory 
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as long as rf <C rj. Only if there are regions in the star where this inequality is reversed somehow could 
damped oscillations occur. 

There may be regions of weak coupling between crustal superfluid and nuclei interspersed among 
regions of strong coupling. If so these regions could, if tied to the strong coupling regions tenously, undergo 
nearly independent oscillations with both long cycle times and insignificant damping. Moreover, there could 
be regions of the core that may undergo long period oscillations if detached from the crust. However, in both 
cases, the effect of nearly independent, slow and persistent oscillations on the portion of the crust where 
superfluid vortex lines are pinned would be minimal, tending to zero in the limit of complete independence. 
Thus, if slow, persistent oscillations can occur somewhere in the star, the chances that one can know about 
them from observations of the rotation rate of that part of the crust where superfluid is pinned strongly 
are remote. In the opposite limit, in which all regions of the star are coupled to one another closely, the 
observed frequencies are averages weighted by moment of inertia, and tend to be dominated by regions of 
high frequency and/or large moment of inertia. 

Under the combined action of external and internal torques, the angular velocity of the crust tends to 
tilt away from alignment with its principal axes. If the external torque is time- independent, or only varies 
on a very long timescale, then ultimately the tilt angle approaches a constant value ~ \N ex \/eI cr Sl, 2 , where 
N ex is the value of the constant external torque, I cr is the typical moment of inertia of the crust, and e is the 
crustal oblateness. If N ex = IQ is the spindown torque, where / is the moment of inertia of the star, then the 
steady state tilt angle is ~ — (I / el cr )il /il 2 , where I is the total moment of inertia of the star. Even though 
— Cl/fl 2 ~ (Cltds) i where t s d is the spindown timescale, is very small (« 5 x 10 _9 P(s)/i s< i(y)) I/eI cr may 
be very large, and 6 could be non-negligible. An amusing side-effect of this tilt is that even an axisymmetric 
neutron star could be a source of gravitational radiation, with an amplitude that can be determined from 
observables (e.g. spindown timescale, period), quantities that can be inferred observationally with varying 
degrees of confidence (e.g. distance) and theoretically determined parameters (e.g. total moment of inertia) 
but does not depend on the oblateness e. Unfortunately, the implied wave amplitudes (strain amplitude 
h < 10~ 30 ) are well below the projected capabilities of the advanced LIGO. 

One key assumption behind this estimate of the steady state tilt angle is contained in the italicized 



word ultimately of the previous paragraph. As was discussed in Section 3.2, the asymptotic value of 9 is 
only attained if the slowest damped mode of the star decays in a time short compared to the timescales on 
which the external torque varies. Practically speaking, this means that if the crustal nuclei and superfluid 
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are closely pinned, then the steady state tilt is approached on the damping time of pd (given by eq. j30[ ) ; if 
this is short compared with the spindown timescale, then the asymptotic value of 6 is reached. Otherwise, 
the tilt could be smaller - somewhere between ~ —(I/<j'I cr )Cl/Vl 2 and —(I/eI cr )il/il 2 - and time variable. 
As was pointed out in Section [}~^, the correct steady state tilt angle, and the evolution toward that angle, 
could not be found from Shaham (1977), where perfect coupling between crust and crustal superfluid was 
assumed. The crux of the solution is in the timescales implied by the imperfection of the pinning. 

We began this paper by proposing to examine perturbations about a particular fixed point of the 
equations governing the rotational dynamics of a neutron star, the one corresponding to equal angular 
velocities of all components lined up along one of the principal axes of the crust (the one with largest 



moment of inertia eigenvalue). We have only wavered from this program briefly, in Section 3.3.1, where 
we considered perturbations about a state in which the rigid crust and crustal superfluid have parallel 
angular velocities with slightly different magnitudes. However, in spite of the constancy of our approach, 
we have uncovered some hints that it may be unrealistic, for when external torques are taken into account, 
the correct fixed points may involve tilts away from principal axes. In a sequel to this paper, we shall 
investigate the implications of time dependent and time independent tilts due to external torques, as well 
as to internal torques we have neglected here. 
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